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We address the problem of finite-size anyons, i.e., composites of charges 
and finite radius magnetic flux tubes. Making perturbative calculations in 
this problem meets certain difficulties reminiscent of those in the problem of 



• pointlike anyons. We show how to circumvent these difficulties for anyons of 



arbitrary spin. The case of spin 1/2 is special because it allows for a direct 



■ application of perturbation theory, while for any other spin, a redefinition of 

the wave function is necessary. We apply the perturbative algorithm to the 
A^-body problem, derive the first-order equation of state and discuss some 
examples. 
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I. INTRODUCTION 



The topology of two-dimensional space allows for existence of anyons [|T],0, particles 
with statistics intermediate between bosonic and fermionic. However, the real world being 
three-dimensional, no particles with inherent anyonic statistics can exist; in reality such 
statistics can only arise effectively by means of an interaction. Remarkably, a composite 
made of a charge e and a magnetic flux $ is an anyon |^, because interchanging two such 
composites multiplies their wave function by exp(z7rQ;), where a = e$/27r, in the spirit of 
the Aharonov-Bohm effect. Such composites arise naturally in Chern-Simons field theory, 
where, the magnetic field being proportional to the charge density, a point charge is at the 
same time a point flux. Quasiparticle excitations in the fractional quantum Hall effect are 
anyons just due to the charge-fiux interaction. 

Since anyons are not fundamental particles but composites or quasiparticle excitations, 
they should have finite size. For example, the size of the FQHE excitations is of the order 
of the magnetic length |Q. Also, in field theory, if the gauge field Lagrangian is the Chern- 
Simons term plus some other term, a point charge generates a magnetic field smeared over 
a finite region. The simplest example is Maxwell-Chern-Simons theory 0,^. Clearly, if 
the radius of the flux tube, which is essentially the inverse photon mass, is much smaller 
than any other distance scale (mean interparticle distance and/or thermal wavelength), 
then the particles are effectively ideal (pointlike) anyons. In the opposite limiting case, one 
has particles in a magnetic field without any change of statistics. Hence there is "distance 
dependent statistics" [0 : If the particles themselves are, say, bosons, they behave like bosons 
when being close together, but like anyons when being far apart. 

One more parameter, the size of the anyons, being present makes quantum mechanical 
problems more complicated. The two-body problem was considered in Refs. , where also 
some general results and conjectures for the A^-body problem can be found. Since even fewer 
results can be derived exactly than for ideal anyons, it is natural to try using perturbation 
theory for small values of the statistical parameter. For ideal anyons, however, perturbation 
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theory gives senseless infinite results in the s-wave sector close to bosonic statistics, due 
to the singular nature of the anyonic interaction. There are several methods to improve 
the situation, which are in fact equivalent to each other ||10| and may be reduced to the 
following: Modify the original Hamiltonian by adding a repulsive contact interaction, 
i.e., consider H'j^ = Hj^ + cYl]'^k^'^i^jk)] this does not affect the exact solutions, but for 



a special value of c all the divergences cancel, leading to the correct finite result pj]-[13 



This value, as was demonstrated in [jT4|, corresponds to assigning to the particles a spin 1/2, 



interacting with the magnetic field inside the flux tube (cf. |T5|). It is also possible (but not 
necessary) to redefine the wave function in such a way that the terms leading to divergences 
cancel already in the Hamiltonian (in addition, also the three-body interaction terms of Hjsi 
then cancel). 

Consider the situation in more detail. Since the nature of the difficulty is connected 
with two-body interaction, it is in fact sufficient to restrict oneself to the two-body problem. 
In the L = sector of the relative Hamiltonian of two ideal anyons there are two linearly 
independent solutions, one regular and the other one singular at r = (r is the interparticle 
distance) but still normalizable. A generic wave function will behave for small r as r"-|-Kr~", 
where k has to be fixed as a boundary condition at r = 0, which corresponds to choosing one 
of the possible self-adjoint extensions of the Hamiltonian [T^. Adding a repulsive contact 



term excludes the singular solution, i.e., fixes k = 0. However, this holds for any value of 
c, while perturbation theory works for one special value only. On the other hand, if the 
singular flux tube is understood as a limiting case of a finite radius one, then normally only 
the regular solution will survive, even with no contact term at all. In fact, the contact term 
has to be regularized itself — and thereby given a physical meaning — and the result (even in 
the singular limit) may depend on the regularization chosen. A simple and natural way is 
to ascribe to the particles a magnetic moment, or spin, which couples to the magnetic field 
inside the flux tube. The solution will then exhibit a continuous dependence on the value 
of the spin, but in the singular limit the problem becomes scale invariant, and therefore the 
dimensionful parameter k, can only tend to oo or to (unless a distance scale is introduced 



by hand [T^, which will not be considered here). That is, the solution can only tend to the 
purely singular one (as we will see, this happens for spin +1/2, where + means attractive, 
i.e., parallel to the flux tube) or to the purely regular one (for any other spin). 

Our present goal is to reproduce these results within (first-order) perturbation theory, in 
order to make possible its application to the A^-body problem. It works best for spin ±1/2, 



considered in detail in Ref. fl^, but does not work directly for any other spin. Analyzing 



the two-body problem, we show how the wave function has to be redefined in order that the 
first-order result be in agreement with the exact solution. We then move on to the A^-body 
problem. Here again, spin ±1/2 is most simple in the sense that it is apparently the only case 
allowing for a cancellation of the three-body terms; however, our algorithm always works in 
the sense that it allows to get rid of the singularities. We apply it to derive the first-order 
perturbative equation of state, using the second-quantized formalism. The limiting cases of 
small and large flux tube radius reproduce the ideal anyon and the mean-field equations of 
state, respectively, and an example illustrates the interpolation between the two. 

II. THE TWO-BODY PROBLEM 

Consider the problem of two finite-size anyons. Let m be their mass, e the charge and 
a the spin (more precisely, the projection of the spin on the z axis, which can be positive 
or negative), a the statistics parameter at infinite distance and e{r) the function describing 
the flux profile (which we assume to be radially symmetric), so that 

$(r-) = ^£(r), Bir) = —e'ir) (1) 
e er 

are the flux through the circle of radius r and the magnetic field at a distance r, respectively; 
one has £:(oo) = 1, and we will assume in the sequel that £:(0) = 0, i.e., there is no singular 



flux at the center (even if the magnetic field still may be singular). As well as in Ref. 
spin does not appear from a relativistic formulation but is introduced by hand. The magnetic 
moment, which couples to the magnetic field, is /z = a— |T^. Then the radial part of the 
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relative Hamiltonian, with a harmonic potential added (the latter, as usually, serves only 

for discretizing the spectrum and is irrelevant for the essence of the matter), is 

^_ 1 / ^2 Id [L-ae(r)]2 2aae'{r)\ 1^^2^2. 
m \ dr"^ r dr r j A ' 

L is the angular momentum, and we assume the "bare" particles to be bosons, whence 

L must be even. We now wish to treat the a dependent terms as a perturbation. The 

difficulties arise for L = 0, when the unperturbed wave functions do not vanish at the 

origin. In fact, it is enough to consider the ground state. For a = 0, its wave function is 

V'o(^) = V rauo exp ^— -mcoT^^ (3) 

and its energy is uj. Assume, for simplicity, that e{r) = 1 for r > R (that is, there is no 
magnetic field outside the circle of radius R), where R is the size of the anyons. Then the 
perturbative correction from the third term of (^ is 

a^e\r) f 1 



mujexp ^— -mcu^r^^ rdr 



= + a UJ - exp — muj r dr 
Jo JR r ^\ 2 J 

— > Q;^ti;(— In i? + finite terms), (4) 

R — ^0 

diverging in the singular limit i? — 0. Therefore, except possibly for some special values 
of a, the ground state of Ti. for small enough R cannot be reproduced perturbatively in a 
straightforward way. (In fact, for R = the exact result, as we will see, is AE = \a\uj for 
any a 7^ +1/2.) 



III. EXACT SOLUTION 

To be able to check the perturbative results, let us first find the exact ground state energy 
of H. Symmetry with respect to a change (a, a) (—a, — o") is present, so it is enough to 
consider a > 0, which is what we will assume, unless otherwise specified. For r > R, the 
wave function with the correct behavior at infinity is 
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where U is the confluent hypergeometric function. The energy E can be determined by 
writing the boundary condition at r = i?, 

ti^ = ^ (a) 
iP{R) R ^ ' 

and substituting first ip->{r) and then ip^{r), the wave function for r < R, thereby obtaining 

two equations connecting E and A. To make the point clear, let us introduce two restrictions: 

(i) a <^ 1, since we are interested in perturbation theory; (ii) q = mujR? /2 <^ 1, since we 

are now concerned about the singular limit, in which g 0. In this approximation]^, the 

first equation in question reads 

E 1 + A - (1 - A)g" 

^ = ' + ^ + A + (l-A)g"- 

Note that substituting the asymptotic form il){r) = + Kr~" into yields the relation 

K = — A)/(l + A). Therefore in the singular limit there are only two possibilities: If 

A ^ — 1 (and 1 + A tends to zero faster than i?^", which will be true for small enough a), then 

the purely singular solution survives and E {l — a)uj, in any other case the regular solution 

survives and E ^ {1 + a)uj. On the other hand, for perturbation theory it is the factor 

in (1^) that is troublesome, because its perturbative expansion = 1 + a In g + " ^ + ■ ■ • 

needs more and more terms as g tends to 0. This factor cancels out for A —>■ ±1, and one 

expects just these two situations to cause no difficulties for perturbation theory. 

To proceed with the exact solution, one has to supply a function e{r) for r < R. The 

simplest choice here is e{r) = r'^/R'^, corresponding to the magnetic field being uniform 

inside the circle of radius R [p|,p|,p!4| . Then (|^) becomes 

1 / rf2 1 ^ \ 1 

H< = - -— + Eo, 8 

m V dr^^ r ar I 4 



^We will, however, continue to refer to the result as "exact", just to stress that it is obtained 
directly from the Schrodinger equation. 
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where 

The corresponding solution — regular at the origin, since there is no singular flux tube — is 



ip^ir) = iFi ( -(1 — ), 1; -mujr j exp \^--mujr ) , (10) 



and substitution into @ yields 



lqE/u + 2aa_ (11) 
a qE/2u}-l ^ ' 



For small enough q one may put A = — 2cr, because q tends to zero faster than Conse- 
quently, the singular solution arises for a = +1/2 (i.e., attractive), the regular one does for 
any other a. Perturbation theory experiences difficulties for any a ^ ±1/2. In particular, 
for the spinless case cr = one gets E/uj = 1 + a(l — g")/(l + g""), i.e., 

/\E = -auj tanh ( - In q) . (12) 



,2 

For g — > 0, this tends to AE = au, but to see this, one has to take into account all orders 



m a. 



IV. PERTURBATION THEORY 



Let us now come back to perturbation theory. Following the idea of |jTT],|T2[, in order to 
get rid of the singularities we redefine the wave function as 



^/,(r) = /(r)V'(r) 



(13) 



where /(r) is to be fixed at our convenience. Then the Hamiltonian acting on ■?/'(r) is 



Ti. = Tio + Til + 7^2, 



(14) 



where 
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n2 



^1 




m 



\ 2 2 

— + -mu r , 
ar I 4 



+ 



a^£^(r) 2aae'{r) 



) 



(17) 



(15) 



(16) 



m 



^2 ^ 



The purpose of the redefinition is to have the dangerous term canceled, by an appropri- 
ate choice of f{r). Asymptotically for r —* oo, when 6{r) = 1, demanding it to be canceled 
by any linear combination of the first two terms of 7^2 will make /(r) be a power function, 
and then each of these two terms will itself be proportional to and therefore should 
be canceled as well. As we will see, the correct way is to demand that all four terms be 
canceled, i.e., that 7^2 = 0, or 



This is the essence of the perturbative algorithm for our problem: Find /(r) from (|T8|), 
choosing the solution which is nonsingular at the origin, and then regard TYi as a perturba- 
tion. This algorithm will yield the correct result in the singular limit. 
Remarkably, for a = =f1/2 there is a universal solution to (0), 



yields the correct first-order result; in the singular limit it is AE = ±au!. In fact, here 
this redefinition is not necessary; acting directly with the original problem (^ is more 
complicated but still possible, because the divergent first-order contribution from the 
term turns out to be canceled by the second-order contribution from the spin term, while 
the first-order contribution from the latter gives the correct answer. 

For arbitrary spin, there is apparently no universal solution. However, use can be made 
of the fact that a -C 1, to construct an approximate one. Coming back to the uniform 
magnetic field model, one has 



r^f'ir) + rf\r) - [aV(r) - 2aare'{r)]f{r) = 0. 



(18) 




(19) 



which is precisely the ansatz considered in |T^. Then one obtains Tii = =F 
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rV"(r) + r/'(r) - - ^]/(r) =0, r < R, ^^^^ 
r^f"{r) + rf'{r) - a^f{r) =0, r > R. 

It is easy to see that /<(r) = 1 — cra^ is the nonsingular solution to the first equation, 
neglecting terms of the order a^. The general solution to the second one is />(r) = Cir" + 
C2r~°, and matching / and df / dr a.t r = R yields 

/W = I (21) 
li-^fe)° + i±^(fr. r>R. 

Now the first-order correction from Tii is 

AE={^o\n,\tljo) 

= 2aaiu- + auj / ) f^e""^ dx. (22 

g ig (1 - 2a)x° + (1 + 2a)g° ^ ^ 

For g -C 1, the first term can be neglected, and in the second term one may replace by 
1 (this is legal unless | lnx| ^ but the main contribution to the integral is given by the 
values X ~ 1), and then the lower limit of integration can be replaced by 0, which yields the 
exact result (|^). 

It is worthwhile to note that the "straightforward" result (^), for the spinless case, is 
correct in a sense; indeed, if in ([l^) one fixes q and goes to the limit a — 0, one does get 
(I). This is natural: There being no singularity, the straightforward perturbation theory 
works for small enough a; but its range of validity shrinks to zero in the limit g — > 0, and to 
get the expression which remains valid in this limit, one has to redefine the wave function 
as described above. 

Concerning the excited states, for the ones with L = and the radial quantum number 
n one has to apply the same algorithm, because its sense is to take care of the short- 
distance behavior of the wave function, which is independent of n. For L 0, perturbation 
theory is directly applicable, because already the unperturbed wave functions vanish at the 
origin, although making the redefinition will do no harm. Thus, the algorithm can be applied 
for all the states, which is essential for the second-quantized formalism. 
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V. THE A^-BODY PROBLEM AND THE PERTURBATIVE EQUATION OF 

STATE 



The main goal of perturbation theory is its apphcation to the A^-body problem. The 
Hamiltonian in our case reads 



N 



HN =—y: 

2m fr{ 



where 



ae^xr 
Mr) = -—^^(0 



(23) 



e 

is the vector potential and B{r) is the magnetic field as in (|T]); rjk = Vj — r^, and Vjk 
The wave function now is to be transformed as 

ipivi, . . . , rjv) = n f\'^3k) V^(ri, . . . , rjv) 

\3<k j 

with /(r) determined as above. Then the Hamiltonian acting on i\} is 



(24) 



(25) 



N 



1 ^ 

2m f^-^ 



92 



d 



f'irjk)rjk d 



+ E 



k^l 



firjk)firji) 



TjkTjl 



TjkTjl 



(26) 



where Eq. (18) has been taken into account. This contains two-body as well as three-body 
interaction terms. If a = ±1/2 then, because of Eq. (p!9|), the three-body terms cancel, 
otherwise they do not. This certainly makes the problem more complicated, however these 
terms are of second order in a and they do not produce any singularities. Therefore in the 
first order one omits them and considers the second and the third terms as perturbation. The 
first-order contribution of the second term, in fact, vanishes, and by virtue of the symmetry 
of the wave function the contribution of the third term can be represented as 

A^(A^- 1) 



AE 



-(^sym|^l(ri,r2)|V^sym), 
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(27) 



where 



^i(ri,r2) 



2 /'(ri2) ri2 d 



(28) 



m firu) ri2 dri 
[cf. (0)]. 

It is now possible to derive the (first-order) perturbative equation of state for finite-size 
anyons near both Bose and Fermi statistics. (In the second well as for bosons outside 

the s-wave sector, the redefinition of the wave function is not necessary but does not change 
the result.) The simplest way is to use the second quantized formalism. The starting point 
is the expression for the first-order correction to the thermodynamic potential 

JO 

(29) 



= ^/^i / ^'^1 ^'^2 7i:i(ri, rs) [{^\r^,f3^)^{r^, A)} {^^{r^, A)V'(r2, A)} 



where |?/'^(ri, /3i)-?/'(r2, /32)} is the one-particle thermal Green function, and the upper/lower 
sign refers to bosons/fermions. Taking into account the explicit form of Tii (note that the 
derivative acts on ip^ only), it is possible to do the spatial integration by parts, and the 
surface term will vanish. There is an expression (see |^ for details) 

oo 

{^\v,, P,)^P{V2, Pi)} = -E(±^)'G(ri,r2;s/3), (30) 
where z is the fugacity and 



s=l 



G(ri,r2;/3) = -^exp 



- r2) 
A2 



(31) 



is the one-particle plane wave thermal Green function in the thermodynamic limit, A 



27r/5/m being the thermal wavelength. Equation (^) then turns into 



^l=«^ECn(±^)" 



(32) 



n=l 



with 



1 



n-l 



1 



I' St 



[°° d ^ 






JO dr 


[ f{r)\ 





1 ± exp 



s + 1 7ir^ 



(33) 
(34) 
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According to Eq. (|T8]), one has 



fir) - 1, 



(35) 



therefore the integral above is in fact proportional to a, and c„ oc 1. Now, one has for the 
pressure 



pp = ~ino + ^i)/v 



A2 



n=l 



(36) 



[Qq = ^{V/X^) Z]^i(='=-2)"/'^^ being the unperturbed thermodynamic potential] and for the 
density 



1 

A2 



=F ln(l :p 2;) - a ^ ncn{±i 



n=l 



neglecting a^, the solution for z is 



2:0 + a(l =F ^0) E riCn{±zoy 



n=l 



where 



zo = ±[1 - exp(TAV)] 



(37) 
(38) 

(39) 

(40) 



is the unperturbed fugacity. Equations ( p6D and (^) explicitly give P as a function of p, 
i.e., the equation of state, which is thus is obtained for any flux profile in terms of integrals. 
Introducing the virial expansion, 

PP = p[l + A,{pX') + AsipX'r + ...], (41) 

the first few virial coefficients are 

^2 = + tt(--ci + C2), 

A3 = 4±a(^ci-2c2 + 2c3), 

17 9 
= a{--ci + -C2 - -C3 + 3C4), 

4 z 

1 1 

= -Trrnrrr ± ^(t^Ci - C2 + Scs - 8C4 + 4C5). 



(42) 



3600 '30 
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VI. EXAMPLES 



Consider some particular cases. We will now drop the condition a > 0. According to 
(p5|), the main term of /(r) at r — oo is proportional to r'"' unless the corresponding C 
vanishes, whence it is proportional to r~'"' [in particular, it is so when = +1/2 and 
e{r) = T]{r)]. In the first case, it is easy to see from (|3^ that 



1±1)^ for i?«A, 



I- St 



\a\ 



for i? > A. 



(43) 



This leads to simple results: 



i? < A : 
i?> A : 



1 ± 1, , 

A2 = +1/4 + 
A2 = +1/4+ 



(44) 
(45) 



and the higher virial coefficients are unaffected. The first equation is the well-known first- 



order result for ideal anyons |20|, the second one is a mean-field result [|21| — indeed, in that 
limit the magnetic field gets smeared over the whole volume. In the second case — that is, 
essentially for attractive spin 1/2 — there is an additional minus sign in the correction terms, 
so that the second virial coefficient can become lower than the bosonic one (cf. |2^). 

The transition between the two limiting regimes can be observed on a simple model where 
an explicit expression for c„ is within reach. Let cr = — 1/2 and e{r) = 1 — exp(— r^/i?^), so 
that the magnetic field B{r) = (2a/ eR^) exp(— r^/i?^) is Gaussian. Then, by virtue of 
rf'{r)/f{r) = e{r), and one gets 

1 



1st 



1 ± 



(46) 



Cn = ± 

n 



2^n^ + 4<^ 



^/n + Vr? + 4<2\ ^f2-n + V^^^ + 4<^ 



^ (n- ^n^ + AnP \ , (2-n - ^ri^ + 4nf^' 
^ +^ ^ 



where 



(47) 
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e = v^f , (48) 

7 is Euler's constant, iIj{x) = r'(x)/r(x). [The same expressions will of course be valid for 
any other a and e{r) that lead to the same /(r), through Eq. ([18|) .] 

At the lowest order in ^ and 1/.^, respectively, the virial coefficients are 

^2 (49) 



5 



3600 



1 1 

36 6e (50) 
1 



16^ 

= -^z:::^:: + a-. 



3600 ■ 60^2 

The Bose and Fermi ^ dependent terms are equal to each other for odd coefficients and have 
opposite signs for even ones, and for ^ <^ 1 they in fact vanish for all odd coefficients but the 
third. Plots of the first five virial coefficients as functions of ^ near Bose statistics, showing 
the intermediate behavior, are displayed in Fig. 1. 

Note that there is no symmetry with respect to a — *■ —a here; indeed, the spin is repulsive 
for a > but attractive for a < 0. 



VII. DISCUSSION 

Let us summarize the main points of our reasoning. In the problem of finite-size anyons 
that was considered, there are no singularities and therefore all results, exact or perturbative, 
are finite. Therefore, had it been possible to calculate the perturbative corrections to all 
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orders (and provided the series converged), the exact result for the problem at hand could 
be obtained without any special treatment. What is achieved by the redefinition of the wave 
function — in a way that actually takes into account the short-distance behavior of the exact 
solution — is the result being obtained at lower order of perturbation theory than it would 
be without this redefinition. For spin 1/2, it is obtained at first order instead of second; for 
any other spin, at first order instead of infinite. 

There is a certain subtlety concerning the transformation of the Hamiltonian. The correct 



singular limit is e{r) = ri{r) (the step function) rather than e{r) = 1 [jl4|, so that the spin 
term in @ becomes a contact term — ^^^5^(r). For a = ^1/2, the function /(r) defined by 
(p!9|) tends in the singular limit to /o('") = r'^°', and it is the equality A In /o(r) = ±27ra;(5^(r) 
that ensures the cancellation of the contact term when the Hamiltonian is redefined |lT0| , p!4| . 
Now, for arbitrary a, the function /(r) defined by (pT]) also tends to r°, but only pointwise; 
an elementary calculation shows that Aln/(r) — > — 47r(ja5^(r), thereby ensuring the correct 
transformation of the Hamiltonian. Again, spin 1/2 is singled out, in the sense that it is in 
this case only that the limiting transition under the operator A is legal. 

The encoding of the short-distance behavior of the two-body wave function in perturba- 
tion theory allows for a perturbative treatment of the A^-body problem. Again, this turns 
out to be most simple for spin 1/2 (because the three-body terms cancel), but still possible, 
in principle, for arbitrary spin. The perturbative equation of state is obtained in an explicit 
form and shows smooth interpolation between the two limiting cases, the ideal anyon one 
and the mean filed one. 

Numerous stimulating discussions with Stephane Ouvry are gratefully acknowledged. I 
am also thankful to Giovanni Amelino-Camelia for useful remarks. Parts of this work were 
done at the theory division of the Institut de Physique Nucleaire, Orsay and at the Centre 
for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, thanks to 
which are due for kind hospitality and support. 
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Figure caption 



Fig. la-d. The virial coefficients 02 through 05 as functions of ^ in the Gaussian model, 
for a = 0.1. 
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Fig. 1a 
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